Semiclassical Analysis of the Wigner 12j Symbol 
with One Small Angular Momentum 
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We derive an asymptotic formula for the Wigner 12 j symbol, in the limit of one small and 11 
large angular momenta. There are two kinds of asymptotic formulas for the 12j symbol with one 
small angular momentum. We present the first kind of formula in this paper. Our derivation 
relies on the techniques developed in the semiclassical analysis of the Wigner 9j symbol [L. Yu and 
R. G. Littlejohn, Phys. Rev. A 83, 052114 (2011)], where we used a gauge-invariant form of the 
multicomponent WKB wave-functions to derive asymptotic formulas for the 9j symbol with small 
and large angular momenta. When applying the same technique to the 12j symbol in this paper, 
we find that the spinor is diagonalized in the direction of an intermediate angular momentum. 
In addition, we find that the geometry of the derived asymptotic formula for the 12j symbol is 
expressed in terms of the vector diagram for a 9j symbol. This illustrates a general geometric 
connection between asymptotic limits of the various 3nj symbols. This work contributes the first 
known asymptotic formula for the 12j symbol to the quantum theory of angular momentum, and 
serves as a basis for finding asymptotic formulas for the Wigner 15j symbol with two small angular 
momenta. 

PACS numbers: 03.65. Sq, 02.30.1k, 03.65.Vf 



I. INTRODUCTION 

This paper derives in full detail one of two kinds of 
asymptotic formulas for the 12 j symbol, in the limit of 
one small and 11 large angular momenta. We will first 
briefly review some of the previous works on the Wigner 
12 j symbol. Its definition and exact formula are de- 
scribed in the textbooks on angular momentum theory 
[H-i]. Although it is used less often than the Wigner 3j- 
, 6j symbols, it has applications in the theory of x-ray 
absorption branching ratios two-photons absorption 
spectroscopy Q, and loop quantum gravity Q. The 12 j 
symbol was first defined in the paper by Jahn and Hope 
Q in 1954. In that paper, they listed two kinds of for- 
mulas for the special values of the 12j symbol when any 
one of its 12 arguments is zero. See Eq. (A8) and Eq. 
(A9) in Q. The complete symmetries of the 12 j symbol 
are given in Q, as a result of the observation that the 
graphical representation of the triangular conditions of 
the 12 j symbol is a Mobius strip. This Mobius strip is 
illustrated in Fig. [1] below. 

We note that Eq. (A8) in |8J corresponds to placing 
the zero argument at the edge of this Mobius strip, and 
Eq. (A9) in [8j corresponds to placing the zero argument 
at the center of the Mobius strip. Thus, we expect that, 
by placing the small angular momentum at the edge and 
center, respectively, of the Mobius strip, we will have two 
kinds of asymptotic formulas for the 12 j symbol with one 
small angular momentum. 

The main theoretical tool we use is a generalization 
of the Born-Oppenheimer approximation, in which the 
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small angular momenta are the fast degrees of freedom 
and the large angular momenta are the slow degrees of 
freedom. The necessary generalization falls under the 
topic of multicomponent WKB theory. See [lfj| - tl2l | for 
the relevant background. The new techniques used in this 
paper are recently developed in the semiclassical analysis 
of the Wigner 9j symbol with small and large quantum 
numbers [13j | . This paper makes extensive use of the 
results from that paper, and assumes a familiarity with 
it. 

In analogy with the setup for the 9j symbol in 13], 
we use exact linear algebra to represent the small angu- 
lar momentum, and use the Schwinger's model to rep- 
resent the large angular momenta. Each wave-function 
consists of a spinor factor and a factor in the form of 
a scalar WKB solution. For the 9j symbol, the scalar 
WKB solutions are represented by Lagrangian manifolds 
associated with a 6j symbol, which have been analyzed 
in 14] to reproduce the Ponzano Regge action [lj|. In 
the problem of the 12 j symbol with one small quantum 
number, the scalar WKB solutions are represented by 
Lagrangian manifolds associated with a 9j symbol. The 
actions for the 9j symbol are presented in [16j. We will 
quote their results in some of the semiclassical analysis 
of the 9j symbol in this paper. 

We now give an outline of this paper. In section [TTl we 
display the spin network of the 12 j symbol in the form 
of a Mobius strip, and decompose it into a scalar prod- 
uct of a bra and a ket. In section [TTTl we define the 12 j 
symbol as a scalar product of two multicomponent wave- 
functions, whose WKB form are derived in section lTVl By 
following the procedure in [l3| , we rewrite the multicom- 
ponent WKB wave-functions into their gauge-invariant 
forms in section fVl In section [VTl we describe the path 
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used in a semiclassical analysis of the Lagrangian man- 
ifolds associated with the 9j symbol by generalizing the 
paths used in [3] ■ We then obtain the action integral as- 
sociated with this path by quoting the results from 
Finally, we calculate the spinor inner products at the in- 
tersections of the Lagrangian manifolds in section I VIII 
Putting the pieces together, we derive an asymptotic for- 
mula for the 12j symbol in section fVIIIl and display plots 
for this formula against exact values of the 12j symbol 
in section IIXI The last section contains comments and 
discussions. 



II. SPIN NETWORK OF THE 12j SYMBOL 

The spin network [4j for the Wigner 12 j symbol 



Jl 32 Jl2 Jl25 
h H J34 J135 
J13 J24 35 36 



(1) 



is illustrated in Fig. [T] In the spin network, each tri- 
angular condition of the 12j symbol is represented by a 
trivalent vertex. The spin network has the shape of a 
Mobius strip. 

J135 




where the small angular momentum is placed at a center 
position at j§. 

One can decompose the spin network of the 12 j symbol 
into two spin network states by cutting 32,33 at the twist, 
and cutting 31,34,35,36 along the center of the strip. Us- 
ing this decomposition, illustrated in Fig. [2j the 12 j sym- 
bol is expressed as a scalar product between a bra and a 
ket, in the Hilbert space represented by the six angular 
momenta ji, ■ ■ ■ ,je- This is explicitly expressed in Eq. 

m. 





FIG. 2. A decomposition of the spin network of the 12j sym- 
bol. 



III. DEFINING THE 12 j SYMBOL 

We use the decomposition of the spin network for the 
12j symbol in Fig. [2] to write it as a scalar product. This 
is equivalent to Eq. (A4) in [8[. We have 



Jl 32 Jl2 Jl25 
33 34 J34 jl35 } = 
313 324 S 5 j 6 



(b\a) 



{[312] [334] b'13] [j'24] [7125] [j'135]} 1 



FIG. 1. The spin network of the Wigner 12j symbol. 



(2) 

where the square bracket notation [•] denotes [c] = 2c+l, 
and I a) and \b) are normalized simultaneous eigenstates of 
lists of operators with certain eigenvalues. We will ignore 
the phase conventions of \a) and \b) for now, since we did 
not use them to derive our formula. In our notation, the 
two states are 



The symmetries of the 12j symbol are associated with 
the symmetries of the Mobius strip, which are given by 
sliding along the Mobius strip and reflecting it about the 
vertical center of Fig. [1] Using these symmetries, any 
position in the center can be moved to any other posi- 
tion in the center, and any position on the edge can be 
moved to any other position on the edge. However, a cen- 
ter position cannot be move to an edge position, or vice 
versa. Thus, there are two inequivalent asymptotic lim- 
its of the 12j symbol with one small angular momentum, 
corresponding to placing the small angular momentum 
at the center or edge of the strip, respectively. In other 
words, we could either place the small angular momen- 
tum at 31,34,35, or jg across the center of the strip, or 
place it at j 2 , 33, 312, J34, J13, 324, jus, or j 135 along the 
edge of the strip. In this paper, we will focus on the case 



\b) = 



h I2 h h S 5 Iq J 12 J34 J125 Jtot 

jl j2 j3 34 S 5 j 6 J12 J34 jl25 



Ii h h h S5 Iq J 13 J24 J135 Jtot 



Jl J2 J3 J4 S5 J6 Jl3 J24 Jl35 







(3) 



(4) 



In the above notation, the large ket lists the operators 
on the top row, and the corresponding quantum numbers 
are listed on the bottom row. The hat is used to distin- 
guish differential operators from their symbols, that is, 
the associated classical functions. 

The states \a) and \b) live in a total Hilbert space of 
six angular momenta Hi ® H2 ® ® V.4 ® He ® Hs- 
Each large angular momentum J r , r = 1, 2, 3, 4, 6, is rep- 
resented by a Schwinger Hilbert space of two harmonic 
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oscillators, namely, H r = L 2 (R 2 ) [13]. The small angu- 
lar momentum S is represented by the usual 2s + 1 di- 
mensional representation of SU(2), that is, H s = C 2s+1 , 
where s — S5. 

Let us now define the lists of operators in Eqs. ([3]) and 
Eq. Q. First we look at the operators I r , r — 1,2, 3, 4, 6, 
J121 ^13! ^24: which act only on the large angu- 
lar momentum spaces T-L r , each of which can be viewed 
as a space of wave-functions ip(x r i,x r 2) for two har- 
monic oscillators of unit frequency and mass. Let a rfl = 
{x ril +ip rii )/^/2, and aJ M = (x r>l -ip rfJ ,)/-/2, ^=1,2, be 
the usual annihilation and creation operators. The op- 
erators I r and J r i are constructed from these differential 
operators a and as follows, 

r = ~2 ^r^r , Jri — ~^ CL r (7idr , (5) 

where i = 1,2,3, and <Ji are the Pauli matrices. The 
quantum numbers j r , r — 1,2,3,4,6 specify the eigen- 
values of both I r and J 2 , to be j r and j r (jr + 1), respec- 
tively. 

The operators J 2 2 , J34, J 2 3 , and J 24 that define in- 
termediate coupling of the large angular momenta are 
defined by partial sums of J r , 

J12 = Ji + J2 , J34 = J3 + J4 • (6) 



Jl 3 =Jl+J 3 , J 24 =J 2 +J4. (7) 

The quantum numbers ji , i = 12, 34, 13, 24 specify the 
eigenvalues of the operators J 2 to be + 1), for i = 
12, 34, 13, 24. See [l7| for more detail on the Schwinger 
model. 

Now we turn our attention to the operator S 2 that acts 
only on the small angular momentum space C 2s+1 . Let S 
be the vector of dimensionless spin operators represented 
by 2s + 1 dimensional matrices that satisfy the SU(2) 
commutation relations 

[Si,Sj] = i €ijk Sk ■ (8) 

The Casimir operator, S 2 = s(s + 1), is proportional to 
the identity operator, so its eigenvalue equation is triv- 
ially satisfied. 

The remaining operators J 2 25 , J135, and Jtot are non- 
diagonal matrices of differential operators. They are de- 
fined in terms of the operators I r , Jri, and Si as follows, 

(JlmW = [J12 + K 2 s(s + 1)]<W + 2J12 ■ S Q(3 , (9) 
(A 2 35 W = [J? 3 + » 2 *(« + l)]5 af3 + 2J 13 • S Q(3 , (10) 

(Jtot)<*/3 = (Jl + J2 + J3 + J4 + 3 6 )S a p + ftSaP- (11) 

These three operators act nontrivially on both the large 
and small angular momentum Hilbert spaces. 



IV. MULTICOMPONENT WAVE-FUNCTIONS 

We follow the approach used in paper [l]| to find 
a gauge-invariant form of the multicomponent wave- 
functions ipa( x ) = (x,a\a) and ipa( x ) — (x,a\b). Let us 
focus on ipa( x )> since the treatment for ip h is analogous. 
We will drop the index a for now. 

Let Di , i = 1, . . . , 12 denote the the operators listed in 
the definition of \a) in Eq. ([3]). We seek a unitary opera- 
tor U , such that Di for alH = 1, . . . , 12 are diagonalized 
when conjugated by U. In other words, 

Ul^D t ) a0 Up u = {k i )^, (12) 

where Aj, i = 1, . . . , 12 is a list of diagonal matrix opera- 
tors. Let (f)^' be the simultaneous eigenfunction for the 
/x th diagonal entries Ai of the operators Ai, i — 1, . . . , 12. 
Then we obtain a simultaneous eigenfunction ipa of the 
original list of operators Di from 

i>^=U afl 4>M. (13) 

Since we are interested in i\) a only to first order in h, all 
we need are the zeroth order Weyl symbol matrix U of U, 
and the first order symbol matrix Ai of Ai . The resulting 
asymptotic form of the wave- function ip(x) is a product 
of a scalar WKB part Be lS and a spinor part r, that is, 

^ (:r) = B(x) e l s{x)/h rW (x, p) . (14) 

Here the action S{x) and the amplitude B{x) are simul- 
taneous solutions to the Hamilton- Jacobi and the trans- 
port equations, respectively, that are associated with the 
Hamiltonians \\ ■ The spinor r M is the p th column of 
the matrix U, 

r^\x,p) = U aM (jc,p), (15) 

where p = dS(x)/dx. 

Now let us apply the above strategy to the 12j sym- 
bol. The Weyl symbols of the operators I r and J r i, 
r = 1,2, 3, 4, 6, are I r — 1/2 and J r i, respectively, where 

Jr ^ ^ ^ Zr^Zrft , Jri — ^ ^ ^ Zr^iip " ) ^LV^rv , (1^) 

and where z r ^ = x rp _ + ip rll and z rM = x r ^ — ip r ^ are the 
symbols of a and a\ respectively. The symbols of the 
remaining operators have the same expressions as Eqs. 
©, 0, ©-CET]), but without the hats. 

Among the operators Di, J 2 25 and the vector of the 
three operators J to t are non-diagonal. By looking at Eq. 
©, the expression for J 2 25 , we see that the zeroth order 
term of the symbol matrix J 2 25 is already proportional 
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to the identity matrix, so the spinor r must be an eigen- 
vector for the first order term J 12 • S. Let r( M )(J 12 ) be 
the eigenvector of the matrix J i2 • S with eigenvalue fxJ\ 2 , 
that is, it satisfies 



(J 12 ' S) Q/3 r^ M) = [iJ 12 T l 



(17) 



where \i = — s, ... , +s. In order to preserve the diagonal 
symbol matrices J\ 2 through the unitary transformation, 
we must choose the spinor t^' to depend only on the 
direction of J12. One possible choice of is the north 
standard gauge, (see Appendix A of [11]), in which the 
spinor 5 a M is rotated along a great circle from the z-axis 
to the direction of J 12. Explicitly, 



only of J12, and since the Poisson brackets {J12, I r } = 
vanish for all r = 1, 2, 3, 4, 6, the second term in Eq. ([22)) 
vanishes. We have 



[A(J r )W = (ir- 1/2) 5p 



(23) 



Similarly, because {J12, Jl 2 } = and {J12, J34} = 0, we 
find 



[A(Ji2)W — ) [-M-^JW — ^34^, 



(24) 



Now we find the symbol matrices A(Ji2s) for the vector 
of operators J125, where 



(3 12) = e^-")*" dW ( 



(18) 



where (#12,012) are the spherical coordinates that spec- 
ify the direction of J12. Note that this is not the only 
choice, since Eq. (|T7l) is invariant under a local U(l) 
gauge transformations. In other words, any other spinor 
t' — e 19 ^ 12 ^ t that is related to r by a U(l) gauge trans- 
formation satisfies Eq. (fT7| . This local gauge freedom is 
parametrized by the vector potential, 



A^=i(r«)t 



83 



12 



(19) 



which transforms as — A^ — Vj 13 (<?) under a 

local gauge transformation. Moreover, the gradient of the 
spinor can be expressed in terms of the vector potential, 
(Eq. (A.22) in as follows, 



[A(Ji25)] M ^ — 

After converting the above operator equation to Weyl 
symbols, we find 



[A(Ji 25 )W (26) 

= Jl2<W - iW^iUa^, J12} + hUa^SapUpv 

= 3 12 S^ - ifirW^Ti"), J12} + nr^*S a pT^ . 

Let us denote the second term above by and use Eq. 
(I2TJ1) . the orthogonality of t, 



00* (v) _ If 
'a 'a — "C" ' 



(27) 



to get 



9 t (m) 



12 



00 

L 12 



J12 x S 

J 12 



r 00 



(20) 



Once we obtain the complete set of spinors t w , fj, = 
— s, . . . , s, we can construct the zeroth order symbol ma- 
trix U of the unitary transformation U from Eq. (|T5|) . 

Now let us show that all the transformed symbol ma- 
trices of the operators in Eq. ([3]), namely, the Aj, are 
diagonal to first order. Let us write A[D] to denote the 
operator U'DU, and write A[D] for its Weyl symbol. 
First, consider the operators I r , r — 1,2,3,4,6, which 
are proportional to the identity matrix. Using the oper- 
ator identity 



[A(/ r )U - ul ll {i r 8 a p)Up v = - u^[u av , i r ] , 

(21) 

we find 



[A(i r )W = ( J r - 1/2)^ - iKUScv I r } , (22) 

where we have used the fact that the symbol of a commu- 
tator is a Poisson bracket. Since U a ^ — is a function 



= -*W*[{rM, Jh} + {tM, J 2i }} 



(28) 



-ih,T^*e kji J : 



(") 



12/v 



<9J 



12 j 



= s(A^ x j 12 ) 4 v + ft^i^ - nr^*s a pA u) , 

where in the third equality, we have used the reduced Lie- 
Poisson bracket (Eq. (30) in [l7| ) to evaluate the Poisson 
bracket {r, Ji} and {t, J 2 }, and in the third equality we 
used dr/Ji = dr/Ji 2 and dr/Ji — dr/Ji 2 from the 
chain rule, and in the fifth equality, we used Eq. (|2"0]) for 
dr/d3\ 2 . Notice the term involving S in T* in Eq. ([28]) 

cancels out the same term in A(Ji2s) in Eq. (|2u| . leaving 
us with a diagonal symbol matrix 



[A(Jl25)]^ = J 



12 



1^2 



■hA$ x J 12 ■ 



(29) 
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Taking the square, we obtain 

[A(Jl26)W = ( J i2 + f-Lh) 2 5^ . (30) 

Finally, let us look at the last three remaining opera- 
tors J tot in Eq. (fTTj) . Since each of the the symbols 
J,, for r = 3, 4, 6 defined in Eq. (fTB")) Poisson com- 
mutes with J12, that is, {J12, Jr} = 0, we find A(J r ) = 
J,. — ihl/Q {Uo (3 12), J r } = Jr, for r = 3,4,6. Using 
A(Ji25) from Eq. (|2"9")l . we obtain 



and r a — with [i — j'125 — ji2- Note that all the 
Hamiltonians except the last three, Jt"ti preserve the 
value of J12 and 3q along their Hamiltonian flows. 

We carry out an analogous analysis for ?p b (x). The 
result is 

^ a (x) = B b (x)e iS "^/ n T b a (x, P ), (34) 

where Sb(x) is the solution to the following 10 Hamilton- 
Jacobi equations: 



[A(Jtot)] //// 
J12 ( 1 



f)+SA^.J l:: 

J12, 



(Ja- 



Je) 



(31) 



Therefore, all Aj, i = 1, . . . , 12, are diagonal. 

The analysis above is completely analogous to those in 
fbij . except that the spinor is diagonalized in the direc- 
tion of the intermediate angular momentum vector J12. 
We see that the procedure in [l3[ generalizes to the case 
of the 12j symbol wave-functions without any complica- 
tion. This is because of the chain rule for differentiation 
and Poisson brackets. See the calculations in Eq. (|28l) . 

Not counting the trivial eigenvalue equation for S 2 , 
we have 11 Hamilton- Jacobi equations associated with 
the Aj for each polarization /1 in the 20 dimensional 
phase space C 10 . It turns out that not all of them are 
functionally independent. In particular, the Hamilton- 
Jacobi equations A(jf 2 ) = J\2^ = C/12 + 1/2)H and 
MJ125) = (J12 + ^h) 2 = (j'125 + l/2) 2 ft 2 are functionally 
dependent. For them to be consistent, we must pick out 
the polarization fi = j'125 — j±2 ■ This reduces the num- 
ber of independent Hamilton- Jacobi equations for S(x) 
from 11 to 10, half of the dimension of the phase space 
C 10 . These ten equations define the Lagrangian manifold 
associated with the action S(x). 

Now let us restore the index a. We express the multi- 
component wave- function ifia(x) in the form of Eq. (|14[) . 



h 


= (h ^ 


-l/2)fi, 


(35) 


h 


= (h H 


-l/2)fi, 




h 


= (is ^ 


-l/2)fc, 




h 


= {k ^ 


-l/2)fc, 




h 


= {je ^ 


-i/2)h, 




I 2 


= (ii3 


+ i/2) 2 h 2 , 




J 2 
J 24 


= U24 


+ i/2) 2 h 2 , 




T (6) 
J tot 


= Jl3 


vh~\ 

1 + — +fiAi 3 xJi3- 
^13. 


h (J 2 + J 4 + Je) = 0. 



Here the spinor r b = rj^ satisfies 



(Jis-S)^ (r^^Ju (r b M )^ 

where z/ = ji 3S - ji 3 . 

The vector potential A13 is defined by 



(36) 



A13 - i{r b ) 



d3 u 



(37) 



Again, note that all the Hamiltonians except the last 
three, j[ b J t , preserve the value of J13 and J 6 along their 
Hamiltonian flows. 



CW = B a {x)e iS ^' n T a a {x,p) . 



(32) 



Here the action S a (x) is the solution to the ten Hamilton- 
Jacobi equations associated with the /x th entries Xf of ten 
of the symbol matrices A", given by 



h = {j x + 1/2)%, 
h = (j2 + l/2)H, 

h = ( j3 + i/2)h, 
h = (jA + i/2)h, 

h = (je + l/2)h, 

JI2 = (ii2 + i/2) 2 h 2 , 



(33) 



J: 



34 
r(«) _ 



34 



'12 



+ l/2) 2 ft 2 , 

LlH 



h A12 x J 12 + (J 3 + J 4 + Je) = . 



V. THE GAUGE-INVARIANT FORM OF THE 
WAVE-FUNCTIONS 

We follow the procedure described by the analysis pre- 
ceding Eq. (69) in 13] to transform the wave-functions 
into their gauge-invariant form. The result is a gauge- 
invariant representation of the wave-function, 



= B a (x)e lS « {x V R [U a (x)r a (x )} 



(38) 



where the action (x) is the integral of p dx starting at 
a point z , which is the lift of a reference point Xq in the 
Lagrangian manifold C®J . The Lagrangian manifold C?J 
is defined by the following equations: 
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h 




-l/2)ft, 


h 


= (h ^ 


-l/2)ft, 


h 


= (j3 "I 


-l/2)ft, 


h 


= (J4 H 


- 1/2)7* , 


h 


= (j'e ^ 


-1/2)71, 


l 2 

J 12 


= (j'l2 


+ l/2) 2 ft 2 , 


T 2 


= (j34 


+ l/2) 2 fi 2 , 


Jtot 


= Jl + 


- J 2 + J 3 + J 4 + Je = 



(39) 



The rotation matrix U a {x) that appears in Eq. (I38[) is 
determined by the SO (3) rotation that transforms the 
shape configuration of J 12 and J6 at the reference point 
z o = {xo,p(xq)) on C® J to the shape configuration of J 12 
and J6 at the point z = (x,p(x)) on L 9 J . Here J12 and 
J6 are functions of z and are defined in Eq. (|16[) . 

Similarly, the multicomponent wave-function for the 
state \b) has the following form, 

i> b {x) = B b (x) e< 3{x ^ h [U b (x) T b (x )} , (40) 

where the action SJ (x) is the integral of p dx starting at 
a point that is the lift of xo onto the Lagrangian mani- 
fold C b . The Lagrangian manifold C b 3 is defined by the 
following equations: 



fcth component. The amplitude 0,^ and the Maslov index 
Hk are the results of doing the stationary phase approx- 
imation of the inner product without the spinors. Each 
rotation matrix U® k is determined by a path -y a (° fe ) that 
goes from z to z% along C 9 J , and U b k is similarly de- 
fined. The formula Eq. (|42[) is independent of the choice 
of Zfc, because any other choice z' k will multiply both 
and U" J by the same additional rotation matrix which 
cancels out in the product {U b k )^U^ k . 

The above analysis is a straightforward application of 
the theoretical result developed in [l]|. We present the 
detail of this analysis to show that the procedure outlined 
in [l3[ does generalize to higher 3nj symbols, such as the 
12 j symbol. 



VI. THE LAGRANGIAN MANIFOLDS AND 
ACTIONS 

We now analyze the Lagrangian manifolds C 9 J and C b J , 
defined by the Hamilton- Jacobi equations Eq. (|39[) and 
Eq. (|4*Tj) . respectively. We focus on C?J first, since the 
treatment for C 9 ^ is analogous. Let k : $5^ — > A 5j denote 
the projection of the large phase space $5^ = (C 2 ) 5 onto 
the angular momentum space A$j = (M 3 ) 5 , through the 
functions J r j, r = 1,2,3,4,6. The first six equations, 
I r = j r + 1/2, r = 1, 2, 3, 4, 6 fix the lengths of the five 
vectors |J r | = J r , r = 1,2,3,4,6. The three equations 
for the total angular momentum, 



h 


= (jl ^ 


-l/2)ft, 


h 


= C?2 ^ 


-1/2)71, 


h 


= (h 4 


-1/2)71, 


h 


= (34 -t 


-1/2)71, 


h 


= (je -t 


-l/2)7i, 


J 13 


= (ii3 


+ l/2) 2 7i 2 , 


T 2 
J 24 


= (j24 


+ l/2) 2 ft 2 , 


Jtot 


= Jl + 


J 2 + J 3 + J 4 + J 6 



(41) 



The rotation matrix U b (x) that appears in Eq. (|4TJj) is 
determined by the 5*0(3) rotation that transform the 
shape configuration of J 13 and J6 at the reference point 
z o = { x q,p(xq)) on C 9 J to the shape configuration of J13 
and J6 at the point z — (x,p(x)) on C b . 

Taking the inner product of the wave-functions, and 
treating the spinors as part of the slowly varying ampli- 
tudes, we find 

(b\a) = e lK Y, ^ exp{i[S 9 J(z k ) - Sf' (z k ) - W 7r/2]/fc} 

k 

(U° k T b (z )) f (U° a k T a ( Zo )). (42) 

In the above formula, the sum is over the components 
of the intersection set Aik between the two Lagrangian 
manifolds C 9 J and C b J . The point is any point in the 



Jtot = Ji + J 2 + J3 + J4 + J 6 = , (43) 

constrains the five vectors Ji, i = 1, . . . , 6 to form a close 
polygon. The remaining two equations 

Jf2 = (ii2 + l/2) 2 fi 2 , (44) 
J34 - (J34 + l/2) 2 7i 2 , (45) 

put the vectors Ji,J2 into a 1-2-12 triangle, and put 
the vectors J3, J4 into a 3-4-34 triangle. Thus, the vec- 
tors form a butterfly shape, illustrated in Fig. [3J This 
shape has two wings (Ji, J2, J12) and (J3, J4, J34) that 
are free to rotate about the J12 and J34 edges, respec- 
tively. Moreover, the Hamilton- Jacobi equations are also 
invariant under an overall rotation of the vectors. Thus 
the projection of C^J onto the angular momentum space 
is diffeomorphic to U(l) 2 x 0(3). 

The orbit of the group U(l) 5 generated by I r , 
r = 1,2,3,4,6 is a five-torus. Thus C^J is a 
five-torus bundle over a sub-manifold described 
by the butterfly configuration in Fig. [3] Alto- 
gether there is a U(l) 7 x SU(2) action on C 9 J . 
If we denote coordinates on U(l) 7 x SU(2) by 
(tpi, ip 2 , tpa, tl>4, 1^6, #12,034,'"), where u e SU(2) and 
where the five angles are the 47r-periodic evolution 
variables corresponding to (I\,l2,l3-,Ii,Ie, J12, J34)! 
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■h 



FIG. 3. The configuration of a point on C 9 J , projected onto 
the angular momentum space A5 j , and viewed in a single R 3 . 



respectively, then the isotropy subgroup is generated by 
three elements, say x = (2tt, 2tt, 2tt, 2tt, 2tt, 0,0,-1), 
y = (0,0,27r,27r,27r,27r,0, -1), and z = 
(27r,27r,0,0,27r,0,27r, -1). The isotropy sub- 

group itself is an Abelian group of eight elements, 
(Z2) 3 = {e, x, y, z, xy, xz, yz, xyz}. Thus the manifold 
C 9 J is topologically U(l) 7 x SU(2)/(Z 2 ) 3 . The analysis 
for C 9 h 3 is the same. 

Now it is easy to find the invariant measure on C 9 J and 
C? h 3 . It is dipx A d^ 2 A d^ 3 A dV>4 A dip 6 A d8 12 A dd 34 A du, 
where du is the Haar measure on SU(2). The volumes 
Va of C 9 J and Vb of C b ° with respect to this measure are 



V A = V B = - (4tt) 7 x 16tt 2 = 2 15 tt 9 , (46) 
8 

where the 1/8 factor compensates for the eight-element 
isotropy subgroup. 

We now examine the intersections of C 9 J and C b 3 in de- 
tail. Because the two lists of Hamilton- Jacobi equations 



(|39f and (|4"Tj) share the common equations I r — j r + 1/2, 
r = 1,2, 3, 4, 6, the intersection in the large phase space 
3>5j is a five-torus fiber bundle over the intersection of the 
projections in the angular momentum space A^j. The 
intersections of the projections in A$j require the five 
vectors J r , r = 1, 2, 3, 4, 6, to satisfy 

|J r | = J r , XX = ' ( 4? ) 

r 

|Ji + J2I = J12 , IJ3 + J4I = J34 , 

|Jl+J 3 | = Jl3, |J 2 + J 4 | = ^24- 

A nice way of constructing the vectors satisfying Eq. (|47|) 
follows the procedure given in the appendix of [18| , which 
was generalized to apply to the symmetric treatment of 
the 9j symbol in llal. For completeness, we summarize 
the construction in [16} using the unsymmetrical labeling 
of the 9j symbol in this paper in the next few paragraphs. 

The construction uses the Gram matrix G of dot prod- 
ucts among the four vectors Ji, i — 1,2,3,4. Some of 
the dot products are given by the length of the vectors 
Ji, i = 1,2,3,4,6, and the intermediate couplings Jj, 
i = 12, 34, 13, 24. In particular, the diagonal elements 
are Jf, i = 1,2,3,4, and some of the off-diagonal ele- 
ments are given by 



Jl 


J 2 = 






- Jl) , 


(48) 


J 3 


J 4 = 


2 (J34 


-4 


- Jl) , 


(49) 


Jl 


J 3 = 




-Jl 


- Ji) , 


(50) 


J 2 


J 4 = 


2 [ 24 


-Jl 


- Jl) ■ 


(51) 



Let us denote the remaining two unknown dot products 
by x = Ji • J4 and y — J2 • J3. We have 



G 



J l 

2 ^ 12 — Jl — J 2 ) 
f( T2 _ 72 _ j2 \ 
2\ J 13 J l J 3) 



\ 



If /2 

2 \ J i2 



2 (Jii 



J 2 
J l 



J 2 



J 2 , 



Jl) 2(^13 



J&) 5(^34 



Jl-Ji) 



J2 

- J 2 



1( j2 _ 72 _ j2 \ 
2\ J 24 J 2 '->!> 



Jl) 



5(^34 



J 3 J I) 
T 2 

J A 



(52) 



The unknown dot products x and y can be solved from a 
system of two equations. The first equation follows from 
Eq. (|43|) . Moving Jq to the other side, and taking the 
square, it becomes 



J 2 - (Ji + J 2 + J 3 + J 4 ) 2 (53) 
= Jl + Jl + Jl + Jl + 2Ji • J 2 + 2J 3 ■ J 4 
+2Ji • J 3 + 2J 2 • J 4 + 2x + 2y , 



which gives us a linear relation between x and y, 



x + y : 



■J4 



Ji- 



■Ji 



Jl 



- V 



- 1 2 - 

J 34 



T 2 - J 2 ) 
J 13 J 24> 



_ (54) 
This is the same equation as Eq. (6) in [16(, except for 
the relabelling of the vectors. The second equation comes 
from the fact that the Gram matrix of the dot products 
between any four vectors in R 3 has a zero determinant. 
That is, 
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P(x,y) = \G\ = 0. 



(55) 



This constitutes a second equation for x and y. Substi- 
tuting the linear relation Eq. (|54|) into Eq. (l55l leads to a 
quartic equation Q(x) — 0, which we can use to solve for 
x. Then we can use Eq. (f5*4"|) to solve for y. In general, 
we find two sets of real solutions of (x, y) — (xi, yi) and 
(x, t/) = (x2,y2). See [16| for more detail. 

For each set of solutions of (x,y), we obtain all the 
dot products among the first four vectors. Assuming all 
the diagonal sub-determinants of order 3 of the Gram 
matrix in Eq. (152[) are positive definite, we can follow the 
procedure outlined in the appendix of [IH to obtain the 
vectors. Let G3 be the first diagonal 3x3 sub-matrix of 
G. We use its singular decomposition to determine the 
vectors Ji, J2, J3. We can then find J 4 from the known 
dot products between 3 i: i = 1,2,3 and J4. Finally, we 
obtain Jq from 



-(Ji+ J 2 + J 3 + J4) 



(56) 



Once we have Ji, i = 1,2,3,4,6, we add them up pair- 
wise to find the intermediate vectors Ji, % = 12, 34, 13, 24. 
This completes the construction of all nine vectors in M 3 . 




FIG. 4. The configuration of a point on the intersection In 
set, projected onto the angular momentum space Asj, and 
viewed in a single R 3 . 



The construction of the vector configuration above not 
only gives explicit solutions for all the vectors at the in- 
tersection of C 9 J and C b J , we also find that there are 
generally two distinct solutions of dot products (x\,yi) 
and (X2, 2/2)- This implies that the solution set consists of 
two sets of vector configurations that are not related by 
an 0(3) symmetry. Thus the solution set of Eq. f|47jl in 
Asj consists of four disconnected subsets, each diffeomor- 
phic to SO(3). These four sets can be grouped into two 
pairs according to the values of the dot products (xi,yi) 
and (x 2 ,y 2 )- 

The intersections in $ 5j - are the lifts of the intersec- 
tions in A5J. Therefore, the intersection of C 9 J consists 



of four disconnected subsets, where each subset is a five- 
torus bundle over SO(3). Let us denote the two sets 
corresponding to (xi,yi) by Iu,Ii2, and denote the two 
sets corresponding to (x2,j/ 2 ) by Ii\,Iti- The vector 
configuration for a typical point in l\\ is illustrated in 
Fig. [4] Each intersection set is an orbit of the group 
U(l) 5 x SU(2), where U(l) 5 represent the phases of the 
five spinors and SU(2) is the diagonal action generated 
by Jtot ■ 

The isotropy subgroup of this group action is Z2, gen- 
erated by the element (27T, 2-7T, 27t, 27r, 2-7T, —1), in coordi- 
nates (tpi, ip2, ip3, "4>4, "067 u ) for the group U(l) 5 x SU(2), 
where u G SU(2). The volume of the intersection man- 
ifold In, Ii2, l2i, or 122, with respect to the measure 
dipi A dip2 A A dip4 A dips A du, is 



1 



13^7 



Vi = 2( 47r ) x 167r = 2 " ;i 



(57) 



where the 1/2 factor compensates for the two element 
isotropy subgroup. 

The amplitude determinant is given in terms of a de- 
terminant of Poisson brackets among distinct Hamiltoni- 
ans between the two lists of Hamilton- Jacobi equations 
in Eqs. (|39|) and (|4T|) . In this case, those are (J12, J31) 
from Eq. fl3g and (J 13 , J 24 ) from Eq. (0TJ). Thus the 
determinant of Poisson brackets is 

\J\2i -hz} {J12, J24} 
{J34, J13} {J34, J24} 



1 



J12J23J13J24 
1 

J12J23J13J24 



Vl23 
V341 

1^123^432 



^214 
V432 



V2I4V34I 



where 



V ijk = J 4 • (Jj x J fe ) . 



(58) 



(59) 



The amplitude ttk in Eq. (|42[) can be inferred from Eq. 
(10) in [16|. In the present case, each fifc has the same 
expression J7. It is 



(27«)Vj VJ12J23J13J: 



24 



VVaVb ^1^123^432 - V214V341I 

(2™)2 13 tt 7 V.A2J23J13J24 



2 15 7T 9 



•\Zl^L23V432 — V214V341 



W J12J23J13J24 



27T\/|Vi23^432 — V214V341I 



(60) 



We now outline the calculation of the relative phase 
between the exponents S a (zi2) — Sb{zi2) and S a (zu) — 
Sb(zn), which can be written as an action integral 

S' (1) = (S a (zi 2 )- S b (z 12 ))- (S a (zu)- Sb(zn)) = <j> pdx 

(61) 
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around a closed loop that goes from zn to Z12 along C 9 J 
and then back along C b ° . 



Hamiltonian flows of (1%, I2, 13, 1 a, Is)- This constitutes 
the last path from p' to p. 




FIG. 5. The loop from a point p G In to q £ I12 along & , 
and then to q' G /12 along 7i2, and then to p' e in along £^ J , 
and finally back to p along In. 



We shall construct the closed loop giving the relative 
phase by following the Hamiltonian flows of various 
observables. This loop consists of four paths, and it is 
illustrated in the large phase space $5^ in Fig. [S] The 
loop projects onto a loop in the angular momentum space 
A$j , which is illustrated in Fig. [6] We take the starting 
point p € In of Fig. [5] to lie in the five-torus fiber above a 
solution of Eq. (|4"T1) . Its vector configuration is illustrated 
in Fig. IDJa). 

First we follow the 3\ 2 how and then the J34 flow to 
trace out a path that takes us along C 9 J from a point 
p in In to a point q in /12. Let the angles of rotations 
be 2^12 and 2034, respectively, where 0i2 is the angle 
between the triangles 1-2-12 and 12-34-6, and 034 is the 
angle between the triangles 3-4-34 and 12-34-6. These 
rotations effectively reflect the triangles 1-2-12 and 3- 
4-34 across the triangle 12-34-6, as illustrated in Figs. 
IHlJa) andjfjb). In addition, the triangle 13-24-6 is also 
reflected across its own plane. Thus, all five vectors J r , 
r = 1, 2, 3, 4, 6, are reflected across the triangle 12-34-6. 

Next, we follow the Hamiltonian flow generated by — jg- 
Jtot along /12 , which generates an overall rotation of all 
the vectors around — j6- Let the angle of rotation be 
2(/) 6 , where 6 is the angle between the triangles 12-34- 
6 and 13-24-6. This brings the triangle 13-24-6 back to 
its original position. However, the triangle 12-34-6 is now 
rotated to the other side of triangle 13-24-6, as illustrated 
in Fig. IDJc). This corresponds to the point q' in Fig. [SJ 

To bring the triangle 12-34-6 back to its original po- 
sition, we follow the J43 flow and 3^ flow along C b ° . 
Let the angle of rotations be 20i3 and 2024, respectively, 
where 0i3 is the angle between the triangle 1-3-13 and 
the triangle 13-24-6, and 024 is the angle between the 
triangle 2-4-24 and the triangle 13-24-6. These rotations 
effectively reflect all the vectors across the triangle 13-24- 
6. Thus we arrive at a point p' € In, where the points p 
and p' have the same projection in the angular momen- 
tum space A5J. This is illustrated in Figs. [Ha) and^d). 
Thus the two points p and p' differ only by the phases of 
the five spinors, which can be restored by following the 




FIG. 6. The loop from Fig. [5]projected onto a loop in As-,-, as 
viewed in a single R 3 . 



To summarize the rotational history in the angular mo- 
mentum space, we have applied the rotations 



fllsQ'lS, 20 13 )i?24(j24> 20 24 )i?(-j 6 , 20 6 ) (62) 
-R34(j34, 20 3 4)-Rl2(jl2, 20i 2 ) , 

where R12 acts only on Ji and J2, -R34 acts only on J 3 
and J4, i?i3 acts only on Ji and J3, i?24 acts only on 
J2 and J4, and R(—je, 20g) acts on all five vectors. The 
corresponding SU(2) rotations, with the same axes and 
angles, take us from point p in Fig. [5] to another point p' 
along the sequence p q —¥ q' — > p' . 

To compute the final five phases required to close the 
loop, we use the Hamilton- Rodrigues formula [l9j . in the 
same way as Eq. (46) in [lj]. Let us start with vector Ji. 
The action of the rotations on this vector can be written 



#(ji3, 20 13 )-R(-j6,20 6 )ii(ji2, 20 12 )Ji = Ji 



(63) 



By inserting an edge Ji6 = Ji + J6 as in part (c) of 
Fig. [3 we split the angle 06 that appears in the mid- 
dle rotation in Eq. (|63|) into two internal dihedral angles 
06a and 06b, of the tetrahedrons in Figs. [TJ^a) and[7jb), 
respectively. 

Then the rotations in Eq. (jBSJ) become 
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FIG. 8. The angles r is the angle between the normals of 
the adjacent triangles sharing the edge J r , where the normals 
are defined by the orientation of the triangles shown. This is 
essentially Fig. 2 in [T^| . with an unsymmetrical labeling of 
the 9j. 



FIG. 7. Decomposition of the angles 4>i an d 06 into dihedral 
angles in two tetrahedra. 



i?(ji 3 , 2<t> 13 )R(-h, 20 6 )fl(ji2, 2^i 2 ) 

= [RQ 13 , 201 3 )iZ(-j6, 206a)] [R(-k, 206fc)i?(jl2, 201 2 )] 

= iZ(ji,20i o )JZ(ji,20i 6 ) 

= i?(ji,20i), (64) 

where we have used the Hamilton-Rodrigues formula 
twice in the second equality. In the third equality, we 
used the fact that 0i — <j>i a + 0i o , where the angles 0i a 
and 0ib are internal dihedral angles for the tetrahedra in 
Figs.[7]Ja) andlT^b), respectively. Thus, we find that the 
product of the three rotations in Eq. (1551) is i?(ji,20i), 
where 0i is the angle between the triangle 1-2-12 and the 
triangle 1-3-13. We can lift the rotation Eq. (|64p up to 
SU(2) with the same axis and angle. Its action on the 
spinor at p is a pure phase. To undo this pure phase, 
we follow the Hamiltonian flow of I\ by an angle — 20i, 
modulo 2-7T. 

Similarly, we can find the rotations acting on J2, J3, J4, 
and J 6 , and proceed to calculate the action integral as in 
[l4T ]. Instead of completing the derivation of the action 
integral using our unsymmetrical labeling of the 9j sym- 
bol, we will quote the result from Eq. (12) in [16|. It is 
given by 

=2J2JMV , (65) 

r 

where ipl 1 ^ = ir — r is the external dihedral angle be- 
tween the normals of the two triangles adjacent to J r . 
The orientations of the triangles are defined in Fig. [8] 
The sum is over r = 1,2,3,4,6,12,34,13,24. The rela- 
tive action integral that corresponds to the other solution 
(x 2 ,y 2 ) of Eq. (gTJ) is 



5^=2^J r 4 2 ', (66) 

r 

which has the same expression as Eq. (|SS^I , but we should 

note that the angles ipr^ are different from ipr , because 
the vector configuration has a different set of dot prod- 
ucts. As in [l(|, we pick to correspond to the root 
in which — n < ip r < n, and pick to correspond to 
the root in which < tp r < tt. 

Altogether, the asymptotic formula for the 97 symbol 
when all j's are large is given by Eq. (1) in [16j, which 
we reproduce here: 

3i h 312 ) 

33 k 334 > (67) 

313 hi 35 J 

1 cos(£W ) 

+ 1 — sin(S< 2 >) . 

It is found from Eq. (17) and Eq. (18) in [16| that, 
when the configuration goes to its time-reversed image, 
that is, when all the vectors reverse their directions, 
the actions transform accor ding to -> -S« and 
-> -S^ + 2tt(J21 =1 >) + 9tt. As a result, the two 
terms cos(S^) and sin^ 2 )) in the 9j formula (|67[) are 
invariant under time-reversal symmetry. In the asymp- 
totic formula (1501) for the 12j symbol that we will derive 
below, the additional phases generated from the spinor 
products will break this time-reversal symmetry. 

Putting the amplitudes ft from Eq. (|50"]) and the rela- 
tive actions and into Eq. we find 
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(b\a) = e iKl - 



V J12J34J13J24 



2ir^/\V (1) V {1) ~V (1) V (1) \ 



+e' 



V J12J23J13J24 



2ir\/\V (2) V (2) - V {2) V (2) \ 

z7r V I v 123 v 432 ^214 ^341 1 



(rH 2u ))t(r«(zu)) + e^ (I, ^-/ 2 )/ s (^r^zn)) + (^ 1 V°( 2ll )) 
{r\z 21 ))\r% Z21 )) + e^ i2) -^/ 2 ^ (ll^ r b (z 21 ))\u^ t« (z 21 ) 



(68) 



where the superscripts (1) and (2) are labels used to dis- 
tinguish the first and the second solutions to Eq. (|4"T|). 
Here we have factored out two arbitrary phases e lKl and 
e iK2 f or ^ w0 p a j rs f stationary phase contributions. 

The rotation matrices Us. , i = 1,2, are determined by 
the paths from zn to Zi 2 along C®J . Similarly the rotation 

matrices Ul f , i — 1,2, are determined by the paths from 
zn to Zi 2 along . See Eq. (76) in [l3[ for a similar, 
but simpler, expression for the case of the 9j symbol. 

VII. SPINOR PRODUCTS 

We choose the vector configurations associated with 
zn to correspond to a particular orientation of the vec- 
tors. We put J12 along the z axis, and put 3q inside the 
xz plane, as illustrated in Fig. [S] Let the inclination and 
azimuth angles (9, <f) denote the direction of the vector 
J 13. From Fig. [9l we see that <j> is the angle between the 
(J12, 3%) plane and the (J12, J13) plane. We denote this 
angle hy <p — 4>i2- The inclination angle 9 is the angle 
between the vectors J12 and J13. 

Z 




X 



FIG. 9. The vector configuration at the point zn in 7n. 

The gauge choices for the spinors at the reference point 
211 are arbitrary, and they only contribute a phase that 



can be absorbed into e tKl . To be concrete, since J12 
points in the z direction, we choose the spinor t°(zu) to 
be the (X th standard eigenvector for S z , that is, 

t£(zii) = <W- (69) 

For the spinor r fc (zn), we choose it to be an eigenvector 
of J13 • S in the north standard gauge, that is, 

T b Jz n ) = e^-^ dl a (9) , (70) 

where (<pi2,0) are the spherical angles of J13 in a refer- 
ence frame where J12 is in the z-direction, and Jq is in 
the xz-plane. See Fig. [5] We denote the azimuthal angle 
by </>i2, because it is also the angle at 3 12 between the 
(J12, J13) plane and the (J12, 3q) plane. 
Taking the spinor inner product, we obtain 

(r b (z n ))t(r a (z n )) = e^-")*. d^(9) . (71) 

To evaluate the other spinor product at Z12, we need 
to find the rotation matrices Ua^ and , which are 
generated from paths j a and 7b from z\\ to z\2 along 
C 9 J and , respectively. 

We choose the path j a to be the path from p to q 
generated by the 3\ 2 flow and the J| 4 flow, which are 
illustrated in Fig. [5] in the large phase space, in Fig. EJa) 
in the angular momentum space. This path contains no 
flow generated by the total angular momentum, so 

C/«=l. (72) 

We choose the path 7^ to be the inverse of the path from 
q back to p along CJ in Fig. \5\ which contains only one 
overall rotation around — jg. Thus 

C/ b (1) =*7(j 6 ,20 6 ). (73) 

The rotation associated with U9~' is illustrated in Fig. 
[SJb). It effectively moves J13 to its mirror image J' 13 
across the 12-34-6 triangle in the rrz-plane, which has 
the direction given by (—(^12, 9). Thus UbT b (zu) is an 
eigenvector of J' 13 • S, and is up to a phase equal to the 
eigenvector of J' 13 • S in the north standard gauge. Thus, 
we have 
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[t4 ( M^ll)]a 



< ( a -"^" d s va (6) , (74) 



where 7?i3 is a holonomy phase factor equal to the area of 
a spherical triangle on a unit sphere; see Fig. [TU] There- 
fore, the spinor product at the intersection I12 is 



{ui x) r\z 11 ))\U a T a {z 11 )) 



//// 



(*)■ 

(75) 




Ti = 



lK1 V J12J3AJ13J2A 



(77) 



X cos 



71" A/IV123V432 — V214V341I 

— + A 1 I 013 I + ^013 

Here H12 is another holonomy for the J12 vector, and 
the angle i3 is the angle between the (J13, J12) 
plane and (J13, J&) plane. Because the quantities 
tyi, 012, 013, -^12, -Hi 3 depend only on the geometry of the 
vector configuration, and are independent of and v, we 
conclude that the argument in the cosine must be linear 
in (i and v. Equating the two arguments of the cosine 
in Eqs. ([7B| and ([77|) . we find that this linear term is 
(M012 + ^013)- Using the Maslov index \i\ = from [l(|, 
we find 




\J J12J34J13J24 



v 123 v 432 



yW V W 
v 214 v 341 



(78) 



■ M0i2 



4%) , 



where we have put back the superscript (1). Through an 
analogous calculation, we find 



FIG. 10. The phase difference between two gauge choices 
can be expressed as an area around a closed loop on the unit 
sphere. 



Let us denote the first term in Eq. (f68|) by T\. Sub- 
stituting the spinor inner products of Eqs. (I7T1) and ((75)) 
into Eq. (|68p . we find that T\ is given by 



e lK2 VJ12J34J13J24 



t/( 2 )t/( 2 ) _ T/WT/V 
^123 r 432 r 214 r 341 I 



r(2) T/ (2) 



x sm 



(V 2 



A 2 ) 

M012 



A 2 ) 

^013 



(79) 



Ti 



e lKl \J J12J34J13J24 



K\/ IU123U432 



X COS 



o(l) W 



V214V341 

M012 



(76) 



#13 
"2" 



312 



Using a different choice of the reference point and 
paths, we can derive an alternative expression for the in- 
ner product, and eliminate the term if 13. Let us choose 
a new reference point Z\\ to correspond to an orientation 
in which J13 is along the z-axis, and 3q lies in the x-z 
plane. We choose the path -f a to go from p to q' along 
the first two paths in Fig. [5j and we choose 7 b to be the 
inverse of the last two paths that goes from q' back to p 
in Fig. [5] Through essentially the same arguments, we 
find 



VIII. ASYMPTOTIC FORMULA FOR THE 12j 
SYMBOL 

From the definition in Eq. ([2]), we see that the fac- 
tor ([ji2][j34][ji3][j24]) 1 ^ 2 in the denominator of Eq. fl2J 
partially cancels out the factor ( J12T34 J13 J24) 1 / 2 from 
Ti and T 2 in Eqs. ([75)1 and (f79"|) . respectively, leaving a 
constant factor of 1/4. Because the 12j symbol is a real 
number, the relative phase between e lKl and e lK2 must 
be ±1. Through numerical experimentation, we found it 
to be +1. We use the limiting case of js = s = from 
Eq. (A9) in Q to determine the overall phase convention. 
This determines most of the overall phase. The rest can 
be fixed through numerical experimentation. Putting the 
pieces together, we obtain a new asymptotic formula for 
the 12j symbol with one small quantum number: 
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Jl 32 312 3125 
33 34 J34 Jl35 
313 324 S j 6 



(-1)" 

47rV(2ji25 + l)(2ji35 + l) 



v 123 v 432 "214 v 341 1 



(IK 

13 I 



(80) 



T/ (2) T/ (2) _ T/ (2) T/ (2), 
^123^432 V 214 V 3H\ 



: sin 



(5( 2 )+M 2 2 ) +M 2 3 ) ) 



As mentioned above, the additional terms from the spinor 
product break the time-reversal symmetry. Thus, it is 
essential that S^ 1 ' and are evaluated at the config- 
urations in which V ~ Jq ■ (J12 x J13) < 0, and not at 
their mirror images. 

Here, the indices on the d-matrix are given by /i = 
J125 — ji2 and v — j'135 — j'13. They are of the same order 
as the small parameter s. The phases S' 1 ' and are 
defined in (p5|) . and the V's are defined by 



Vijfe — Ji • (Jj x J/j) . 



(81) 



The angles 0x2 and 0i3 are internal dihedral angles 
at the edge J12 and J13, respectively, of a tetrahedron 
formed by the six vectors J12, J13, J24, J34, J6, and J2'3, 
where J2'3 — J3 — ^2- This tetrahedron is illustrated in 
Fig. [IT] The angle 9 is the angle between the vectors J12 
and J13. The explicit expression for the angles </>i2, 4>\s, 
and 9 are given by the following equations 



!>12 = 7T — cos 



^13 



cos 



cos 



(J12 X J13) ■ (J12 X J 6 ) 

|Jl2 X J i3 | I J12 X J 6 | 
(Jl 3 X J12) ■ (J13 X J 6 ) 
|Jl3 X J 12 1 |Jl3 X J 6 | 
Jl2 ' Jl3 
J12J13 



(82) 
(83) 
(84) 



IX. PLOTS 

We illustrate the accuracy of the approximation Eq. 
(|5D|) by plotting it against the exact 12 j symbol in the 
classically allowed region for the following values of the 
j's: 



31 32 312 3125 
33 34 334 3135 
jl3 324 S 5 j 6 

The result is shown in Fig. | 



51/2 59/2 21 22 
55/2 53/2 27 26 
27 25 1 j 6 

From the error plot in 



Fig. n"47 a). we see that the agreement is excellent, even 
for these relatively small values of the j's. 
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FIG. 12. Comparison of the exact 12j symbol (vertical sticks 
and dots) and the asymptotic formula (|80p in the classically 
allowed region away from the caustics, for the values of j's 
shown in Eq. (|85jl . 



Since the asymptotic formula (|80|) should become more 
accurate as the values of the j's get larger, we plot the 
formula against the exact 12 j symbol for another exam- 
ple, 



FIG. 11. The angles 0i2 and (^13 are internal di- 
hedral angles in the tetrahedron with the six lengths 
Je, J12, J34, J13, J24, J2'3, where J 2 '3 = J3 — J2- The angle 6 
is the angle between J12 and J13. 



211/2 219/2 91 92 
= { 205/2 223/2 107 108 
99 93 2 j 6 

(86) 

in the classically allowed region away from the caustic in 
Fig. H31 These values of the j's are roughly four times 
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FIG. 13. Comparison of the exact 12j symbol (vertical sticks 
and dots) and the asymptotic formula (|80[) in the classically 
allowed region away from the caustics, for the values of j's 
shown in Eq. 
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FIG. 14. Absolute value of the error of the asymptotic formula 
(|80[1 for (a) the case shown in Eq. (|85[) , and (b) the case shown 
in Eq. I|86[) . The error is defined as the difference between the 
approximate value and the exact value. 



those in Eq. (|85l) . The errors for this case are displayed 
in Fig. Il"47 b). By comparing Figs. Il~47 a) and ll~4T b) . we 



can conclude that the error scales with the j's. 



X. CONCLUSIONS 

In this paper, we have derived an asymptotic formula 
of the 12j symbol with one small angular momentum, 
generalizing the special formula of the 12 j symbol, Eq. 
(A9) in 8]. By looking at the other special formula for 
the 12j symbol, Eq. (A8) in we can guess that the 
other asymptotic limit of the 12 j symbol will involve the 
semiclassical analysis of the trivial 9j symbol, which re- 
duces to a product of two 6j symbols. We will present 
that result in a future paper. 

The analysis of the 12 j symbol in this paper is a nat- 
ural extension of the analysis of the 9j symbol in [l3j]. 
Based on the calculations in these two papers, we can 
summarize our steps in finding asymptotic formulas for 
the 3nj symbols with small and large quantum numbers. 
First, we ignore the small quantum numbers and any of 
the large quantum numbers that involve the indices of the 
small ones. For instance, in this paper, j$ — S5 is small, 
so we ignore j'5, j'125, and j'135. The remaining relevant 
large quantum numbers determine the Lagrangian man- 
ifolds. Once we fix the Lagrangian manifolds, the scalar 
WKB parts of the wave-functions can be derived from 
a semiclassical analysis of these Lagrangian manifolds, 
following the procedure in [3, E3] • The spinor parts of 
the wave-functions at the intersection points of the La- 
grangian manifolds arc determined by the path used to 
calculate the action integral in the semiclassical analysis. 
Finally, taking the inner product of both the scalar part 
and the spinor part of the wave-functions, we can derive 
an asymptotic formula for the 3nj symbol with small and 
large angular momenta. 

In general, we note that the asymptotic limits of a 3nj 
symbol with one small angular momentum is expressed 
in terms of the geometry associated with the asymptotic 
limits of a 3mj-symbol, where m = n — 1. Since the 
Wigner 15j symbol is used extensively in loop quantum 
gravity and topological quantum field theory, we suspect 
that there are deeper, and more geometrical interpreta- 
tions of these approximate relations of the 3nj symbol in 
their various semiclassical limits. 
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